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Chapter  1 


INTRODUCTION 

The  inception  of  cavitation  is  caused  by  a  reduction  of  local 
pressure  that  leads  to  the  formation  of  bubbles  in  a  liquid.  The 
occurrence  of  cavitation  is  usually  inevitable  on  underwater  bodies 
moving  at  high  speeds,  and  designers  are  constantly  challenged  to  limit 
and  suppress  its  inception.  Two-dimensional  slots  along  a  wall  are,  at 
times,  unavoidable  in  the  designing  process.  Low  pressure  regions 
inside  the  slot  can  lead  to  cavitation.  Specifically,  cavitation  usually 
occurs  in  the  center  of  the  vortex  found  inside  a  two-dimensional 
rectangular  slot.  The  fluid  mechanics  of  flow  in  a  two-dimensional  slot 
must  be  understood  before  further  progress  can  be  made  toward 
predicting  the  inception  of  cavitation  that  occurs  within  the  slot. 

The  purpose  of  the  research  presented  here  is  to  study  the  fluid 
mechanics  of  flow  over  a  rectangular  slot  on  a  flat  plate  and  to  document 
the  magnitude  and  location  of  the  internal  vortex  with  changes  in  slot 
geometry.  The  parameters  of  critical  importance  to  the  flow  include  the 
ratio  of  boundary  layer  thickness  (based  on  99%  of  the  freestream 
velocity)  and  slot  width  (8/B),  the  slot  span/width  (S/B),  and  the 
depth/width  (D/B)  ratios.  Of  interest  were  cases  where  the  boundary 
layer  was  thicker  than  the  slot  width  (i.e. ,  8/B  was  greater  than  unity,'.  In 
order  to  study  the  fluid  mechanics  inside  the  slot  without  cavitation 
actually  occurring,  air  rather  than  water  was  c.ncsen  as  the  fluid  mecium. 


During  the  tests  the  ratio  of  boundary  layer  thickness  (8)  to  slot  width  was 
fixed,  while  S/B  and  D/B  (from  now  on  referred  to  as  the  aspect  ratio, 
AR=D/B)  were  varied. 

1 . 1  Previous  Work 

Much  of  the  previous  research  on  flow  in  two-dimensional  slots 
was  found  to  be  inapplicable  to  the  present  work.  However,  the  problem 
of  slot  cavitation  seemed  to  be  the  incentive  for  the  first  research  done. 
Bennett  (1949)  studied  the  applied  problem  concerning  the  downstream 
effects  of  a  clearance  slot  around  the  top  hatch  of  the  water  tunnel  test 
section  at  the  Ordnance  Research  Laboratory  at  The  Pennsylvania  State 
University.  He  generalized  the  problem  to  include  the  fluid  mechanics 
past  a  slot  used  to  guide  the  gates  of  dams.  His  research  on  slot  fluid 
mechanics  was  performed  in  a  wind  tunnel.  Boundary  layer  thickness, 
freestream  velocity,  model  width,  model  depth  and  the  downstream 
corner  radius  were  the  vanables  in  his  study.  Bennett's  emphasis  was 
the  fluid  mechanics  of  the  slot  entrance  and  the  downstream  corner,  not 
the  flow  field  contained  within  the  slot.  This  follows  since  his  concern 
was  to  suppress  cavitation  at  the  downstream  corner  which  would  disturb 
the  flew  within  the  water  tunnel  test  section.  It  seems  that  the  cavitation 
that  prcbaDiy  occurred  inside  the  slot  was  net  of  concern  or  not  yet 
cosenyed  since  one  of  his  conclusions  was  that  "depth  had  veiy  little 
effect"  (Bennett,  p.  65,  1949)  on  the  minimum  pressure  found  at  the 


ccwnstream  corner. 
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A  large  amount  of  the  research  performed  on  the  two-dimensional 
slot  was  aimed  at  the  slot's  effect  on  the  passing  flow.  This  includes  the 
effects  of  lew  subsonic  oscillations  due  to  shallow  slots  (Sarohia  1377). 
Sarohia  seemed  particularly  concerned  with  the  effect  of  the  shear  layer 
across  the  slot  entrance  on  the  resulting  oscillations.  Similar  work  was 
performed  in  a  water  channel  by  Rockwell  and  Knisely  (1979)  but  they 
investigated  the  effect  of  the  downstream  eege  on  the  unsteadyness  of 
the  flow.  Both  of  these  investigations  had  a  laminar  boundary  layer  at  the 
upstream  corner  of  their  respective  slots. 

Military  applications  found  a  need  for  research  on  high  subsonic 
and  low  supersonic  flow  across  shallow  rectangular  slots.  Clark  and 
Kauffman  (1980)  presented  research  on  this  type  of  flow  for  Mach 
numbers  of  0.6  to  3.0.  Their  emphasis  was  on  the  aerodynamic  and 
aero-acoustic  effects  due  to  a  slot.  However,  they  conclude  that  more 
work  would  be  required  to  estimate  the  presure  distributions  in  the  slot 
(Clark  and  Kauffman  1980). 

Several  investigations  have  examined  the  three-dimensionality  of 
the  flow  inside  a  rectangular  slot.  Mauil  and  East  (1963)  stated  that  the 
flow  in  a  rectangular  slot  was  not  two-dimensional  anywhere  inside. 

They  showed  that  a  regular  three-dimensional  flow  developed  in  the 
rectangular  slot  .  The  regular  three-dimensional  type  of  flow  is 
reenforced  by  the  research  of  Hiwoda,  Mabuchi  and  Masaya  (1985). 
However,  their  description  of  this  flow  is  quite  different  and  more  complex 
than  that  of  Mauil  and  East.  This  could  be  due,  in  part,  to  the  fact  that 
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Maull  and  East  had  a  laminar  boundary  layer  at  the  upstream  corner  of 
tneir  slot  while  Hiwada,  Mabuchi  and  Masaya  had  a  small  (less  than  their 
slot  width)  turbulent  boundary  layer  at  the  same  location.  It  is  worth 
noting  that  Hiwoda,  Mabuchi  and  Masaya  concluded  that  "the  three- 
dimensional  flow  behavior  is  independent  of  the  span  length"  (p.  89, 
1985).  This  does  not  infer  that  the  flow  is  two-dimensional;  rather,  their 
emphasis  is  on  a  cell  type  structure  that  repeats  down  the  model  span. 

A  mainstay  of  most  of  the  above  work  was  that  performed  by 
Roshko  (1955).  His  wind  tunnel  experiments  showed  the  unsteady 
nature  of  the  vortex  inside  of  the  two-dimensional  slot.  He  addressed 
many  of  the  issues  of  the  two-dimensional  slot  including  the  internal 
vortices,  induced  drag  and  transport  of  mass  into  the  slot.  However,  his 
boundary  layer  thickness,  located  1.5  inches  (38.1  mm)  upstream  of  the 
slot  entrance,  was  measured  to  be  a  third  of  the  slot  width. 

The  experiments  most  pertinent  to  the  present  research  came  from 
Ward  (1973),  and  Holl  and  Billet  (1986).  Both  investigations  were 
performed  with  turbulent  boundary  layer  thicknesses  equal  to  or  greater 
than  the  slot  width.  Holl  and  Billet  performed  laser-Doppler  velocimeter 
tests  in  a  wind  tunnel  and  obtained  velocity  profiles  of  the  internal  slot 
flew  field.  A  minimum  pressure  coefficient  of  the  slot  was  obtained  by 
integrating  the  measured  velocity  profile  assuming  a  classical  vortex 
structure.  This  pressure  coefficient  (Cp)  defined  as 


yielded  a  very  low  minimum  pressure  coefficient.  This  is  in  contrast  to 
their  internal  slot  vortex  cavitation  index  data  from  water  tunnel  tests, 
which  gave  very  high  indices  of  cavitation;  where  cavitation  index  is 
defined  as 
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If  the  cavition  observed  during  the  water  tunnel  test  was  purely 
vaporous  then  the  cavitation  index  data  would  be  equal  to  the 
analogous  negative  pressure  coefficient  values  obtained  from  wind 
tunnel  tests.  Ward  also  documented  cavitation  index  values  of  the  same 
order  as  Holl  and  Billet.  Both  also  noted  that  peak  cavitation  occurs  at  an 
aspect  ratio  of  approximately  1 .0.  The  large  difference  between  the 
water  tunnel  cavitation  index  data  and  the  wind  tunnel  negative  pressure 
coefficients  led  Holl  and  Billet  to  question  whether  the  observed 
cavitation  is  purely  vaporous  or  partially  gaseous.  Ultimately,  Holl  and 
Billet  concluded  that  "the  determination  of  an  incipient  cavitation  scaling 
relationship  for  slots  remains  an  unsolved  problem  because  it  is  very 
difficult  to  separate  the  effects  of  gaseous  cavitation  inception  from  the 
minimum  pressure  coefficent.  .  ."  (Holl  and  Billet,  p.  11,  1 986). 


The  primary  objective  of  the  present  research  was  the  direct 
pressure  measurement  of  the  pressure  field  inside  a  two-dimensional 
slot,  with  a  large  turbulent  boundary  layer  at  the  slot  entrance.  All 
measurements  were  to  be  conducted  with  a  non-cavitating  slot,  for  the 
purpose  of  investigating  the  fluid  mechanics  associated  with  cavitation 
inception  in  slots.  Toward  this  end,  the  following  specific  objectives  are 
listed: 

1.  Find  a  repeatable  and  consistent  means  of  generating  a 
turbulent  boundary  layer  along  the  wall  upstream  of  the  slot 
such  that  its  thickness  is  larger  than  the  slot  width. 

2.  Determine  the  effect  of  model  span  on  the  two-dimensionality 
of  the  flow  field  by  means  of  flow  visualization. 

3.  Collect,  by  means  of  flow  visualization,  qualitative  information 
on  the  internal  fluid  mechanics  of  a  rectangular  slot  and  the 
mass  transfer  across  the  slot  entrance  as  a  function  of  slot 
aspect  ratio. 

4.  Collect  location  and  strength  data  on  the  internal  slot  vortices 
as  a  function  of  aspect  ratio  by  means  of  direct  pressure 
measurement. 

5.  Investigate  the  reliability  of  the  direct  pressure  measurements 
by  means  of  flow  visualization.  Namely,  determine  if  the 
pressure  probe  utilized  for  direct,  pressure  measurement  alters 
the  internal  slot  fluid  mechanics. 


Chapter  2 


FACILITIES  AND  EQUIPMENT 

The  experiments  reported  here  were  conducted  at  The 
Pennsylvania  State  University  Mechanical  Engineering  Subsonic  Wind 
Tunnel  Facility.  The  tests  were  performed  on  a  suction-type  wind  tunnel 
with  a  test  section  of  dimensions  96  inches  (2.44  m)  in  length,  38  inch 
(0.97  m)  width,  and  1 2  inches  (0.31  m)  in  height.  The  top  of  the  wind 
tunnel  was  partitioned  in  several  8  inch  (0.20  m)  sections  spanning  the 
tunnel  width  which  could  be  lifted  from  the  top.  It  was  in  this  way  that 
access  was  gained  to  the  inside  of  the  test  section  and  also  how  models 
were  mounted  (see  Figure  1).  The  tunnel  had  a  speed  range  of  zero  to 
50  ft/s  (15.2  m/s).  Wind  tunnel  speed  was  controlled  by  a  frequency 
controller  (Parametrics,  Model  6075)  feeding  power  to  the  wind  tunnel 
motor  which  also  shielded  the  motor  from  external  electrical  interference. 

The  equipment  used  in  the  lab  includes  a  differential  pressure 
transducer  (MKS  Baratron,  Model  398HD)  and  a  pressure  signal 
conditioner  (MKS,  Model  270B)  with  a  combined  accuracy  of  +  0.08%  of 
the  display  and  an  available  resolution  of  10'6  mm  Hg.  This  unit  had  a 
wide  enough  range  for  use  in  all  pressure  measurements  including  the 
dynamic  head  and  local  pressures  in  the  slot.  A  x-y  traverse  system 
(Maxwell  Electronics,  Model  SMC-102A)  with  a  senai  port  for  computer 
control  was  used  for  the  boundary  layer  hot  wire  placement.  This 
traversing  equipment  offered  positioning  accuracy  down 


to  1/4000  of  an  inch  (6.3  pm).  The  hot  wire  implemented  a  common 
Wheatstone  bridge  circuit  and  a  low  pass  filter  (Sathianathan  1986). 

The  hot  wire  sensing  area  size  was  0.00015  inch  (4  pm)  in  diameter  and 
0.050  inch  (1.27  mm)  long  (TSI,  Model  1261A-T1.5).  A  personal 
computer  (AT&T,  Model  6300)  with  hard  disk  was  used.  This  computer 
was  equipped  with  a  I/O  parallel  port  used  for  the  traverse  control  and  a 
16  channel  12  bit  data  acquisition  board  (Metra  Byte,  Model  Dashl6F). 
The  data  acquisition  board  was  configured  to  measure  voltages  from 
zero  to  five  volts;  any  signals  below  or  above  this  range  are  cut  off  and 
interpreted  as  zero  or  five  volts,  respectively.  The  software 
accompanying  the  data  acquisition  board  and  its  configuration  offered 
some  restrictions.  For  example,  only  12  bits  were  available  for  analog  to 
digital  conversion,  which  means  the  measured  voltages  were  stored  in 
integer  form  from  0  to  4095.  This  limited  the  resolution  of  the  stored  data 
and  in  the  0  to  5  volt  configuration  this  translated  to  1 .22  millivolt 
resolution. 

The  smoke  wire  technique  (Corke  et  al.  1977)  was  employed  for 
visualization  of  the  flow.  A  0.005  inch  (0.13  mm)  diameter  nichrome  wire 
was  used.  The  equipment  included  a  50  volt  DC  power  supply  (Heathkit, 
Model  IP-20)  and  strobotach  (General  Radio  Company,  Model  1531).  A 
customized  timing/synchronization  controller  unit,  built  for  the  wind  tunnel 
laboratory,  allowed  for  precise  millisecond  timing  and  control  of  the 
smoke  and  strobe  (Sathianathan  1985).  The  smoke  was  photographed 


with  a  35  mm  camera  (Nikon  with  50  mm  lens)  using  black  and  white, 
400  speed  film  (Kodak  Tri-X  pan,  24  exposure,  400). 


:icw  Visualization  Model 


The  slot  models  were  constructed  out  of  half-inch  plexiglass  to 
allow  for  easy  construction  and  transparency  sc  that  photos  could  be 
taken  of  the  two-dimensional  plane  of  interest.  The  internal  width  of  the 
model  was  set  by  the  boundary  layer  measurement  at  1.5  inches  (38.1 
mm).  This  investigation  includes  aspect  ratios  up  to  three,  so  the  model 
was  made  accordingly  deep  enough,  and  this  can  be  seen  in  the  sketch 
in  Figure  2.  To  maintain  easy  control  of  the  aspect  ratio  a  movable  base 
was  used  in  all  of  the  models  and  could  be  adjusted  from  the  outside  by 
means  of  two  1 /4-inch  threaded  rods. 


Only  one  model  was  used  in  the  flow  visualization.  However,  the 
model  was  later  modified  with  two  dividers  to  take  t^e  reduced  span 
photographs.  These  dividers  were  set  to  reduce  the  span  from  38  inches 
(0.97  m)  down  to  8  inches  (0.20  m)  centered  around  the  plane  of  interest. 

To  prevent  glare  and  increase  the  quality  of  the  pictures,  the  inside 
of  the  model  was  covered  with  flat  black  paper  on  four  of  the  five  sides. 
The  fifth  side  was  keot  clear  for  the  camera.  To  get  light  into  the  two- 
dimensional  plane  wnere  smcne  was  placed,  a  half-inch  portion  of  the 
pacer  was  removed  from  the  upstream  and  downstream  sides  of  the 
mccel.  This  haif-mch  slice  can  be  seen  in  the  photograph  of  the  model 


Based  on  ihe  results  of  the  flow  visualization,  the  span  of  the 
original  model  could  be  reduced,  without  significant  affect  on  the  flow,  to 
8  inches  (0.20  m)  centered  around  the  two-dimensional  plane  of  interest. 
A  new  model  was  then  constructed  for  pressure  measurements.  Half¬ 
inch  plexiglass  was  again  used  as  the  model  material  and  the  pressure 
measurement  model  had  the  same  internal  dimensions  as  in  Figure  2. 
The  differences  in  this  model  can  best  be  seen  in  the  sketch  and 
photograph  of  the  pressure  model  found  in  Figures  4  and  5.  To  pass  the 
pressure  probe  down  the  span  the  two  end  plates  were  machined  with  a 
slot  coming  up  from  the  bottom.  Two  manual  x-y  traverses  (Daedal, 

Model  4526)  were  mounted  to  each  face  as  seen  in  the  photograph.  The 
two  traverses  allowed  for  placement  of  the  pressure  probe  perpendicular 
to  the  two-dimensional  plane  and  offered  movement  in  the  x-y  directions. 
Prior  to  mounting  in  the  wind  tunnel  the  probe  was  manually  centered  in 
the  field  through  which  the  traverses  could  travel.  This  was  done  with  a 
pair  of  verniers  and  the  absolute  dimensions  of  the  probe  relative  to  the 
origin  in  Figure  2  were  recorded  for  future  use. 

The  pressure  probe  was  made  from  stainless  steal  hypodermic 
needle  tubing  with  a  1/1 6  inch  (1 .6  mm)  O.D.  and  0.009  inch  (0.23  mm) 
wail  thickness.  The  shorter  model  span  was  necessary  to  prevent 
bcwmg  of  the  probe  in  the  middle  of  the  span.  The  probe  diameter 
represents  4,2%  of  the  model  width.  The  probe  could  net  be  moved  any 
closer  than  0.3  inches  (7,6  mm)  from  the  slot  entrance  and  interference  cf 


the  model  base  plate  and  the  traverse  field  occurred  for  aspect  ratios  less 
than  0.5.  As  a  result  no  measurements  were  possible  below  AR=0.5. 

The  probe  had  four  holes  drilled  through  the  walls  and  distributed  evenly 
around  the  circumference  near  the  model  midspan.  The  holes  were 
placed  on  the  probe  in  this  fashion  so  that  the  local  static  pressure  could 
be  measured.  It  was  assumed  that  the  average  of  the  four  pressure  taps 
would  approximate  the  local  static  pressure. 

The  model  was  sealed  at  its  seams  with  tape;  however,  the  groove 
through  which  the  pressure  probe  passed  could  not  be  sealed  as  easily. 
A  common  plastic  wrap  was  used  for  this  purpose.  Rubber  bands  were 
used  around  the  tube  itself  while  a  thin  stiff  plastic  plate  was  screwed 
down  to  the  end  plates  to  pinch  the  plastic  wrap  against  the  model.  By 
supplying  enough  slack  in  between  the  probe  and  pinch  points  the  probe 
was  free  to  move  without  restriction  in  the  traverse  field  while  maintaining 
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PROCEDURE  AND  RESULTS 
3.1  Boundary  Laver 

A  thick  boundary  layer  was  desired  along  the  wind  tunnel  wall 
upstream  of  the  slot.  Thus,  the  first  task  involved  in  the  project  was  to  find 
a  means  of  tripping  the  boundary  layer  to  a  consistent  thickness  (5).  The 
most  flexible  means  of  boundary  layer  control  was  found  to  be  rough 
grain  sandpaper.  Specifically,  our  purposes  were  served  by  32  inches 
(0.81  m)  of  number  12  grain  sandpaper.  That  is,  the  first  32  inches 
(0.81  m)of  the  wind  tunnel  test  section  were  covered  with  sandpaper  to 
trip  the  boundary  layer.  During  the  boundary  layer  measurement  more  or 
less  sandpaper  could  be  added  or  removed  in  8  inch  (0.20  m)  sections 
from  40  inches  (1.02  m)  down  to  8  inches  (0.20  m).  This  is  illustrated  in 
Figure  1 .  Note  that  for  these  boundary  layer  surveys  the  slot  was  not 
installed  in  the  wind  tunnel. 

The  boundary  layer  was  measured  by  means  of  hot  wire 
anemometry.  The  procedure  was  entirely  computer  controlled  and  the 
software  which  was  written  exclusively  for  boundary  layer  measurement 
can  be  found  in  Appendix  D.  The  procedure  involved  a  calibration  of  the 
hot  wire,  measurement  of  the  profile  velocities,  then  a  second  calibration. 
For  the  calibration,  the  free  stream  velocity  was  measured  by  means  of  a 
pitot-static  probe.  The  hot  wire  and  pressure  measurements  were 
conditioned  to  analog  voltages  which  were  measured  and  recorded  by 
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computer.  To  avoid  aliasing,  the  hot  wire  voltage  signal  was  passed 
through  a  low  pass  filter  set  at  50  Hz.  During  data  acquisition  the 
computer  would  sample  at  1 00  Hz  for  1 0  seconds  for  each  point.  This 
was  dene  for  both  calibrations  as  well  as  for  the  profile  measurements. 

The  computer  was  configured  to  measure  voltages  from  zero  to 
five  volts  and  stored  the  data  in  a  binary  file.  The  data  was  later 
converted  to  real  numbers  with  each  record  containing  1000  readings. 

The  hot  wire  voltage  and  velocity  pairs  from  the  calibration  runs  and  the 
hot  wire  voltages  at  each  position  in  the  boundary  layer  were  averaged. 

The  calibration  files  were  then  used  to  find  a  fifth  order  least-square  fit  of 
velocity  versus  voltage.  The  two  calibration  curves  were  time  weighted 
when  the  profile  voltages  were  matched  with  these  curves.  This  was 
done  to  avoid  problems  caused  by  room  temperature  changes  over  the 
period  of  time  required  for  each  test.  A  typical  profile  test  lasted  15 
minutes,  where  the  temperature  could  change  by  as  much  as  1°  F 
(0.6°  C).  This  software  can  also  be  found  in  Appendix  D. 

For  each  run,  temperature  and  pressure  were  recorded  by  the 
operator.  These  measurements  were  recorded  at  the  beginning  of  each 
of  the  two  calibration  runs.  This  information  was  also  used  for  later 
ncndimensionalization  of  the  boundary  layer  profiles  (e.g.,  to  obtain  a 
value  for  viscosity). 

The  goal  here  was  to  meet  the  given  criterion  set  in  the  project 
proposal  (i.e.  the  boundary  layer  thickness-to-model  width  ratio  (5/B)  was 
to  be  greater  than  unity).  The  boundary  layer  thickness  based  on  99%  of 


the  freestream  velocity  to  be  1.8  inches  (45.7  mm),  measured  along  the 
top  of  the  test  section  at  the  location  of  the  upstream  corner  of  the  slot  (i.e. 
65  inches  (1 .65  m)  downstream  of  the  tunnel  entrance)  with  32  inches 
(0.81  m)  of  number  12  sandpaper  at  the  front  and  a  freestream  velocity  of 
15.75  ft's  (4.8  m/s).  The  model  width  was  then  set  at  1.5  inches  (38.1 
mm)  which  yielded  6/B  =  1 .2.  For  the  remainder  of  the  tests  this  same 
configuration  of  sandpaper  was  installed,  and  the  boundary  layer 
thickness  just  upstream  of  the  slot  was  therefore  fixed. 

The  measured  boundary  layer  profile  can  be  seen  in 
nondimensionalized  form  in  Figure  6.  The  procedure  used  to  get  the 
wall-friction  velocity  v*  was  as  follows:  First  the  momentum  thickness  (9) 
was  calculated  using  the  definition 


This  was  done  numerically  and  the  computer  code  can  be  found  in 
Appendix  D.  Next  the  Reynolds  number  based  on  momentum  thickness 
was  calculated 


Boundary  Layer  Profile  (32  in.  of  sandpaper,  15.75  fl/s). 


The  skin  friction  coefficient  was  then  calculated  according  to  Kays  and 
Crawford  (1 980)  as 
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=  .0125  (Re)-0'25. 


Finally  the  velocity  scale  was  then  calculated  to  be 

i 
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The  boundary  layer  was  then  plotted  in  traditional  fashion  as  seen  in 
Figure  6.  Note  that  in  this  case  the  positive  y  direction  is  downward  from 
the  top  inside  of  the  wind  tunnel  test  section  into  the  freestream. 

The  solid  live  on  the  plot  was  obtained  from  Spalding's  law-of-the- 
wall  equation, 
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where  k=0.4  and  B=5.5  (White  1974).  This  curve  was  included  for 
comparison  and  is  said  to  be  one  of  the  simplest  and  most 
comprehensive  equations  of  turbulent  Doundary  layers. 

Comparison  of  the  experimental  data  with  the  Spalding  profile  shows  a 
right-hand  shift  of  the  experimental  profile.  Overall,  the  experimentally 
determined  profile  is  a  typical  response  of  a  turbulent  boundary  layer  to  a 


wall  disturbance  (Clauser  1956).  However,  the  reason  for  the  shift  is  stiil 
of  some  concern.  The  most  likely  reason  is  a  lack  of  recovery  time  for  the 
flow  from  the  and  of  the  sandpaper  until  the  mcdellocation.  The 
obvious  correction  for  this  discrepancy  would  be  to  move  the  model 
farther  downstream  in  the  wind  tunnel  test  section.  However,  the  wind 
tunnel  freestream  is  not  very  uniform  in  the  Iasi  third  of  the  test  section,  so 
the  model  could  not  be  moved  downstream  without  sacrificing  flow 
quality. 


3.2  F'ow  Visualization 

The  flow  visualization  portion  of  the  project  had  many  different 
goals.  The  first  was  to  collect  some  qualitative  information  on  the  flow 
structure  inside  of  the  cavity  as  a  function  of  the  aspect  ratio  AR  =  D/B. 
The  second  was  to  determine  if  the  two-dimensional  slot  flow  structure 
was  altered  by  shortening  the  span  from  the  full  model  span  of  38  inches 
(0.97  m)  to  some  shorter  value  that  would  be  more  conducive  to  pressure 
measurement  inside  the  slot.  The  third  goal  was  to  see  how  the  primary 
vortex  center  moved  as  a  function  of  the  aspect  ratio.  The  primary  vortex 
was  defined  as  the  strongest  and  closest  vortex  to  the  slot  entrance. 
These  visualizations  also  served  as  a  guide  for  the  design  of  the 
additional  mocel  needed  for  the  pressure  measurements.  Finally  an 
investigation  of  the  slot  entrance  region  was  to  be  carried  out  with  the 
smoke  wire  flew  visualization  technique.  Namely,  the  goal  was  to 
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investigate  the  mechanism  that  maintains  the  vertex  in  the  cavity  and  to 
see  the  extent  of  mass  transfer,  if  any,  across  the  slot  entrance. 


The  method  used  to  visualize  the  flow  is  called  the  smoke  wire 


method  (Corke  et  al.  1977).  The  source  of  smoke  comes  from  a  thin 
mc.orome  wire  passed  down  through  the  center  of  the  slot.  This  wire  is 
weighted  at  the  bottom  and  kept  taut  with  a  spring  from  the  top.  To 
induce  smoke  in  the  cavity  a  mineral  oil  is  dripped  down  the  wire  leaving 
small  droplets  behind  on  the  wire.  A  high  DC  voltage  is  then  placed 
across  the  wire  which  heats  up  and  burns  the  oil.  The  smoke  created  by 
this  process  is  then  carried  with  the  flow  and  the  resulting  set  of 
streaklines  can  be  followed. 


Creating  the  smoke  in  the  cavity  was  not  difficult.  The  difficulties 
came  in  trying  to  synchronize  the  smoke  and  strobe  flash  in  order  to 
photograDh  the  process.  Fortunately,  a  timing/synchronization  controller 
was  available  which  provided  millisecond  control  of  the  smoke  wire  and 
the  strobe.  The  strobe  was  set  up  at  the  half-inch  ports  on  the  sides  of 
the  mode1  to  provide  an  instantaneous  flash  of  light  after  some  desired 
time  delay.  The  procedure  was:  coat  the  wire  with  oil,  turn  out  the  room 
lights,  open  the  camera  shutter,  switch  on  the  timing  controller,  then  close 
the  shutter  immediately  after  the  strobe  flashed.  The  timing/  synchronize  - 
tion  controller  would  be  preset  to  delay  the  strobe  flash  so  that  the  smoke 
had  time  to  be  drawn  into  the  vortex.  The  delay  time  between  charging 


the  smoke  wire  and  triggering  the  flash  was  varied  to  provide  a  variety  of 
smoke  wire  pictures.  The  above  procedure  was  followed  for  aspect 
ratios  of  0.5,  0.75,  1.0,  1.25,  1.5,  1.75,  2.0,  2.5  and  3.0.  One  entire  roil  of 
film  was  taken  for  each  aspect  ratio  in  order  to  examine  the 
unsteadyness  of  the  flew. 

As  will  be  discussed  later,  the  pressure  measurements  inside  the 
slot  required  a  slot  of  significantly  shorter  span  than  the  flow  visualization 
model  which  spanned  the  entire  38  inches  (0.97  m)  of  the  wind  tunnel 
test  section.  Thus,  it  was  necessary  to  determine  if  a  reduced  span 
would  affect  the  two-dimensionality  of  the  flow.  It  was  reasoned  that  if 
changes  in  the  two-dimensionality  of  the  flow  take  place  with  changes  in 
the  span  this  would  be  most  visible  in  the  deepest  aspect  ratio  case  (i.e. 
AR=3.0).  The  flow  was  first  photographed  with  a  span  of  38  inches 
(0.97  m)  which  when  nondimensionalized  by  model  width  is  equivalent 
to  a  span  aspect  ratio  of  25.3.  Next  the  flow  visualization  was  repeated 
for  a  reduced  span  of  8.0  inches  (0.20  m),  equivalent  to  a  5.33  span 
aspect  ratio.  When  compared,  no  difference  was  found  in  the  primary 
vortex  location,  see,  for  example,  Figures  7  and  8.  Based  on  the  y- 
coordinate  of  the  primary  vortex  center  relative  to  an  origin  at  the 
upstream  entrance  of  the  slot  (see  Figure  3),  the  38-inch  (0.97  m)  span 
model  had  an  average  of  0.99  inches  (25.2  mm)  for  the  primary  vortex 
depth;  while  the  8-inch  (0.20  m)  model  offered  an  average  vortex  center 
depth  of  1.02  inches  (25.9  mm).  The  difference  is  only  about  3%  which  is 
surprisingly  small  considering  the  large  amount  of  scatter  due  to 


unsteadyness  in  both  sets  of  smoke  wire  pictures.  These  measurements 
come  from  the  average  of  approximately  20  pictures  for  each  set  and 
both  sets  had  a  scatter  of  approximately  ±  0.28  inch  (7.1  mm).  An 
additional  error  of  +  0.1  inch  (2.5  mm)  from  measuring,  by  hand,  the 
vortex  center  location  from  the  photographs  makes  the  total  possible 
error  of  vertex  vertical  location  ±  0.38  inch  (9.6  mm).  This  scatter  was 
charactenstic  of  all  of  the  photographs  and  is  most  likely  due  to  the 
unsteadyness  of  the  flew.  Since  the  span  change  had  no  apparent  effect 
on  the  flew  the  rest  of  the  photographs  were  taken  with  the  8-inch  (0.20 
m)  model  for  convenience.  The  pressure  measurement  model  was  then 
designed  for  a  span  of  8  inches  (0.20  m)  as  well. 

The  slot  flew  field  was  visualized  for  the  aspect  ratios  previously 
mentioned.  The  photographs  were  taken  using  the  aforementioned  flow 
visualization  procedure  and  can  be  found  in  Appendix  A  for  all  of  the 
aspect  ratios.  The  pictures  in  Appendix  A  were  chosen  based  on  their 
clarity  and  for  their  ability  to  show  the  flow  structure  with  a  well  defined 
primary  vortex  center.  When  examining  the  photos  it  is  advisable  to  keep 
some  of  the  inherent  problems  of  the  smoke  wire  method  in  mind. 
Specifically,  low  velocities  found  at  the  higher  aspect  ratios  did  not  allow 
for  distinct  streaklines.  This  may  be  due  to  the  diffusion  process  of  the 
smoke  or  possibly  three-dimensional  effects  on  the  smene  in  the  cavity. 
Note  that  as  the  smoke  moved  out  of  the  piane  of  interest  the  smoke  also 
moved  out  of  the  focused  area  of  the  camera.  This  can  be  seen  in  the 
higher  asoect  ratios  of  2.0  and  3.0  where  the  smoke  appears  to  have 


23 


folded  over  itself.  These  higher  aspect  ratios  with  inherently  lower  cavity 
velocities  required  a  larger  delay  between  the  closing  of  the  smoke  wire 
curcuit  and  the  strobe  flash.  This  brings  the  additional  possibility  that  the 
relatively  hot  smoke  may  have  nsen  in  the  cavity  due  to  buoyancy  effects. 
These  effects  could  lead  to  misrepresentation  of  the  vortex  center  and 
cistcrticn  of  the  smoke  image  from  the  actual  vortex  structure. 

Fortunately,  this  is  net  observed  in  the  photographs.  If  buoyancy  of  the 
smoke  was  a  factor,  the  photographs  with  longer  strobe  delays  would 
have  indicated  progressively  higher  vertex  centers;  but  this  was  not 
observed. 

This  method  adequately  showed  the  overall  flow  structure  in  the 
slots.  One  notable  observation  from  the  photographs  was  that  the  higher 
aspect  ratios  had  more  than  one  vortex  in  the  slot.  The  secondary 
vortices  were  found  deep  in  the  slot  and  were  weaker  than  the  primary 
vortices.  This  was  seen  by  the  lack  of  smoke  in  the  secondary  vortex 
relative  to  the  primary.  The  indication  was  that  the  velocities  were  lower 
deeper  into  the  slot.  The  low  aspect  ratios  (less  than  2.0)  did  not  develop 
the  secondary  vortex  but  the  intermediate  aspect  ratios  of  1 .75  and  1 .5 
seemed  to  have  a  relatively  calm  region  deep  in  the  slot. 

The  quantitative  result  data  obtained  from  these  photos  is  the 
approximate  location  of  the  primary  vortex  center  as  a  function  of  aspect 
ratio.  It  is  at  the  vortex  center  that  the  lowest  pressures  are  to  be  found  in 
the  cavity,  and  that  is  where  cavitation  would  most  likely  have  occurred  if 
these  test  had  been  performed  in  water  rather  than  air.  The  smoke  wire 
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photographs  were  viewed  on  a  1:1  scale  with  the  model  so  that  the 
vertex  centers  could  be  measured  easily.  They  showed  that  the  vortex 
center  relative  to  an  origin  at  the  upstream  entrance  oi  the  mode!  (see 
Figure  2)  primarily  moves  in  the  y-direction  with  changes  in  the  aspect 
ratio.  From  these  photos  it  was  too  difficult  to  tell  if  the  vortex  moves 
upstream  or  downstream  in  the  x-coordinate.  The  vertical  center  location 
was  plotted  in  two  nendimensionalized  forms  (by  width  and  by  depth)  in 
Figures  9  and  10.  The  first  figure  with  y  noncimensionalized  by  width 
also  showed  the  error  involved.  This  error  included  the  approximation  of 
±  0.1  inch  (2.54  mm)  resulting  'tom  measurements  off  of  the  frequently 
fuzzy  photographs.  The  majority  of  the  error  came  from  the  scatter  found 
in  the  twenty  or  more  pictures  measured  for  each  aspect  ratio.  The  point 
plotted  is  the  averaged  value  found  for  each  case  and  the  error  bars 
indicate  the  range  of  the  vertical  vortex  position.  Despite  the  growing 
scatter  with  increasing  aspect  ratio,  the  y-location  increases  steadily  up 
to  AR  of  about  2.0  and  then  levels  off.  It  is  significant  that  the  secondary 
vortex  also  first  appeared  deep  in  the  slot  at  an  aspect  ratio  of  about  2.0. 

3.2.2  External  Cavity  Visualization 

The  exact  mechanism  of  sustaining  the  vortex  in  the  slot  is  not 
completely  understood.  Specifically,  is  the  vortex  maintained  solely  by 
the  shear  from  the  flow  passing  over  the  entrance  or  is  there  a  more 
involved  interaction  with  mass  being  transferred  through  the  slot 
entrance?  This  question  was  addressed  by  moving  the  smoke  wire  out 
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of  the  slot  and  placing  it  just  in  front  of  the  entrance  while  pnctcgrachmg 
the  slot  interior.  The  procedure  for  taking  the  photographs  was  the  same 
as  for  the  internal  flew  visualization.  Four  different  ascect  ratios  were 
chosen,  namely  AR  =  1.0,1 .5,  2.0  and  3.0.  Examples  of  the  flow  at  eacn 
asoect  ratio  are  found  in  Appencix  3.  The  pictures  presented  were 
chosen  to  demonstrate  the  contrast  between  the  times  when  smoke  was 
and  was  net  seen  in  the  slot.  An  overall  examination  of  the  photegrapns 
shows  that  about  three  out  of  four  pictures  hac  smoke  in  the  slot.  This 
would  indicate  that  some  mass  transfer  is  indeed  taking  place,  and  that 
shear  at  the  slot  entrance  is  net  the  sole  cause  of  vortex  motion.  The 
other  interesting  observation  is  that  smoke  was  never  seen  near  the 
center  of  the  primary  vortex. 

3.3  Pressure  Measurement 

The  objective  in  this  stage  of  the  project  was  to  measure  the 
pressure  field  in  the  two-dimensional  plane,  midspan  of  the  slot.  In 
particular,  the  location  and  strength  of  the  vortex  center  as  a  function  of 
aspect  ratio  was  desired.  This  task  required  measuring  the  static 
pressure  throughout  the  plane  without  disturbing  the  vortex  structure. 
Again  we  note  that  the  primary  vortex  is  mainly  of  interest  since  it  was 
shown  by  flow  visualization  to  be  the  strongest  vortex  whenever  more 
tKan  one  vortex  is  present.  Direct  pressure  measurement  was 
accomplished  by  passing  a  thin  tube  with  holes  in  its  circumference 


down  the  third  dimension  in  the  model  span.  The  details  of  the  model 
itself  can  be  found  in  Chapter  2. 

The  procedure  used  to  survey  the  pressure  field  in  the  slot  was  as 
fellows:  First  the  pressure  probe  was  centered  in  the  traverse  field  and 
the  x  and  y  coordinates  relative  to  the  absolute  origin  (upstream  corner  of 
the  slot)  were  recorded.  The  relative  x  and  y  coordinates  read  from  both 
traverses  were  also  recorded.  The  software  used  during  the  procedure 
required  the  three  pairs  of  coordinates  for  user  interaction  and  storage  of 
the  data.  This  software  can  be  found  in  Appendix  D. 

Just  prior  to  the  pressure  field  surveys,  the  wind  tunnel  freestream 

velocity  was  measured  by  means  of  a  pitot-static  probe  placed  in  the 
freestream.  The  dynamic  head  (1/2  p  uM2)  was  later  used  to 

nondimensionalize  the  measured  pressure  difference  between  the 
freestream  static  pressure  (PJ  and  the  static  pressure  inside  the  slot  (Pc). 

This  nondimensionalization  provided  a  negative  local  pressure 
coefficient  defined  as: 


This  form  was  convenient  since  it  is  anaiogus  to  a  cavitation  index. 

Following  the  dynamic  head  measurement  the  positive  side  of  the 
pressure  transducer  was  connected  to  the  freestream  static  pressure  tap, 
and  the  negative  side  was  connected  to  the  pressure  probe  inside  the 


slot.  In  this  manner  the  pressure  coefficient  ccuid  be  measured  as 
defined  above.  The  probe  was  traveled  through  a  predetermined 
number  of  points  in  a  square  grid  contained  within  the  traverse  limits. 

Each  data  point  was  stored  in  a  binary  file  along  with  the  correspcncing 
absolute  x  and  y  coordinates. 

The  software  was  written  such  that  either  a  7  x  7  point  grid  with 
1/16  inch  (1.59  mm)  spacing  or  a  14  x  14  point  grid  with  1/32  inch  (0.79 
mm)  spacing  could  be  chosen.  Obviously  the  14x14  point  grid  would 
allow  for  higher  resolution  plots;  however,  one  test  run  with  a  14  x  14 
point  grid  would  take  about  5  hours  compared  to  only  1  1/2  hours  for  the 
7x7  grid.  The  14x14  grid  was  used  for  the  aspect  ratios  of  1 .0, 1 .5,  2.0 
and  3.0,  while  the  7  x7  grid  was  used  for  all  the  others  (AR  =  0.5,  0.6, 
0.75,  0.9,  1.25,  1.75,  2.5). 

The  long  periods  of  time  involved  in  each  procedure  became  a 
concern  in  the  preliminary  tests.  Of  particular  concern  was  the  fact  that 
changes  in  room  temperature  could  effect  the  pressure  measurement 
equipment.  Therefore,  the  pressure  signal  conditioner  had  to  be 
rezerced  periodically  during  each  test.  Specifically,  the  equipment  was 
zeroed  prior  to  starting  a  pass  in  each  row  in  the  grid.  In  other  words,  the 
equipment  was  zeroed  every  seventh  point  in  a  7  x  7  grid  and  every 
fourteenth  point  in  a  14  x  14  point  grid.  The  freestream  velocity  was 
readjusted  after  every  seventh  or  fourteenth  point  as  well. 

Despite  this  extra  effort  temperature  changes  in  the  lab  still  may 
have  had  an  effect  on  the  data.  The  largest  observed  deviations  from  the 
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equipment  during  these  tests  were  on  the  order  of  ±  10‘5  mmHg.  This 
represents  approximately  6%  of  the  minimum  pressure  measurements 
(which  were  around  0.CC018  mmHg). 

Early  measurements  discovered  some  interesting  flow 
characteristics  that  required  procedure  adjustments  from  those  used  in 
the  boundary  layer  measurements.  Specifically,  the  very  unsteady 
nature  of  the  vortex  in  the  slot  required  longer  sampling  times.  In  the  final 
test  runs,  data  were  sampled  at  100  Hz  for  30  seconds.  Ten  second 
periods  were  initially  attempted  but  produced  inconsistent  results. 

Secondly,  to  accurately  measure  these  very  low  pressure  differences,  the 
range  selection  dial  on  the  pressure  signal  conditioner  had  to  be 
adjusted  for  maximum  sensitivity.  This  adjustment  also  resulted  in  more 
efficient  use  of  the  A/D  converter  voltage  range.  In  particular,  for  the 
lowest  aspect  ratios(0.5  and  0.6)  the  0.1  x  range  selection  was  used, 
while  the  0.01  x  setting  was  used  for  all  other  aspect  ratios. 

The  pressures  for  each  point  in  the  grid  of  each  test  were 
averaged  and  the  mean  values  were  nondimensionalized  by  the 
dynamic  head.  This  produced  a  square  plot  of  negative  pressure 
coefficients  in  x-y  coordinates  which  were  nondimensionalized  by  width. 
Plots  of  this  kind  demonstrated  that  the  pressures  were  lower  near  the 
vortex  center.  In  order  to  bring  out  this  information  more  clearly  a 
computer  code  was  written  to  check  and  interpolate  for  specific  pressure 
coefficient  values  along  the  rows  and  also  the  columns  of  the  grid.  The 
computer  codes  used  to  average,  nondimensionalize,  sort,  and 
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interpolate  the  data  can  be  found  in  Appendix  D.  This  procedure 
identified  isobars  in  the  grid  and  the  plots  can  be  found  in  Appendix  C. 

Due  to  the  wide  range  of  pressure  coefficient  values  most  plots 
representing  different  aspect  ratios  had  to  be  treated  individually.  The 
deeper  aspect  ratio  tests  could  be  sorted  along  the  same  set  of  isobars 
but  the  tests  with  AR  less  than  1.5  were  processed  separately.  To  avoid 
confusion  each  plot  contains  its  own  legend. 

Quantitative  information  about  the  location  and  relative  magnitude 
of  the  primary  vortex  can  more  easily  be  seen  with  the  compiled  data 
presented  in  several  new  forms.  For  example,  the  horizontal  location  of 
the  vortex  center  as  a  function  of  aspect  ratio  can  be  seen  in  Figure  1 1 . 

Note  the  lack  of  movement  of  the  primary  vortex  in  the  x-direction  with 
aspect  ratio.  This  supports  earlier  observations  from  the  flow 
visualization.  However,  the  vortex  at  the  aspect  ratio  of  0.5  did  not  stay  in 
the  same  horizontal  location  as  did  the  higher  ones  but  shifted 
downstream.  Additional  tests  were  done  with  the  aspect  ratios  of  0.6  and 
0.9  to  supplement  the  data  in  this  region  and  they  show  the  same  trend. 

It  should  be  noted  that  the  lower  aspect  ratios  do  push  the  limits  of  our 
model  the  most.  The  fact  that  the  pressure  probe  is  taking  up  a  larger 
percentage  of  space  in  the  cavity  does  little  for  the  confidence  of  these 
data. 

The  next  two  plots  on  the  vortex  location  are  similar  to  the  ones 
seen  earlier  from  the  flow  visualization  data.  Figures  1 2  and  1 3  show  the 
change  in  the  vertical  vortex  location  with  aspect  ratio.  Eoth  plots  have 
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the  flow  visualization  data  for  comparison  and  Figure  12  contains  error 
bars  as  seen  on  the  visualization  plot  earlier.  The  pressure  data  error 
was  calculated  based  on  mode!  tolerance  and  repeatability  of  the  results. 
This  error  was  ±  0.1  inch  (2.5  mm)  for  all  aspect  ratios. 

The  pressure  and  visualization  data  are  close  and  consistent 
within  the  error  limits.  The  data  at  aspect  ratio  of  0.5  indicate  that  the 
pressure  probe  never  went  low  enough  to  find  the  true  vortex  center.  For 
the  larger  aspect  ratios,  the  visualization  data  are  consistently  higher 
than  the  pressure  data.  This  may  be  due  in  part  to  buoyancy  of  the 
smoke.  The  effect  of  the  probe  on  the  flow  is  also  a  possible  cause  of 
error.  To  test  this,  a  dummy  pressure  probe  (i.e.  a  1/1 6  inch  (1 .59  mm) 
O.D.  rod)  was  placed  across  the  entire  span  at  various  locations  inside 
the  flow  visualization  model.  Photographs  with  this  rod  in  the  model 
were  taken  using  the  smoke  wire  method.  Aspect  ratios  of  0.5,  1.0  and 
2.0  were  chosen  and  these  photographs  can  be  seen  in  Figures  14-18. 
These  all  show  that  the  flow  is  not  affected  by  the  added  rod. 

Finally  the  strength  of  the  primary  vortex  is  shown  in  Figure  19  as 
a  function  of  aspect  ratio.  This  figure  shows  that  the  pressure  coefficient 
decays  with  increasing  aspect  ratio.  Obviously  this  curve  must  originate 
at  the  origin  so  it  would  seem  that  the  tests  did  not  extend  to  low  enough 
ascect  ratios  to  include  both  sides  of  the  expected  peak.  Unfortunately, 
physical  limitations  of  the  model  precluded  measurements  below  an 
ascect  ratio  of  0.5.  The  overall  error  with  this  pressure  data  is  ±  0.003  in 
terms  of  pressure  coefficient.  This  was  determined  from  the  error  created 
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Chaster  4 

DISCUSSION  AND  CONCLUSIONS 

The  fluic  means". cs  associated  with  the  flew  of  a  two-dimensional 
slot  in  a  hat  cate  "as  ceen  s'udiec:  soecifical'y,  the  case  where  a 
turDu.ent  ocur.cary  ayer  m  ckness  of  t  .2  tiroes  the  slot  width  has  been 
studied.  The  flow  fleic  ~  dscar  of  an  8-.ncn  (0.20  mj  span  slot  was 
visualized  by  means  of  tue  s^cke  wire  method  and  the  flow  field  was 
subjected  to  cirect  treasure  measurement.  The  flow  visualization 
qualitatively  cescnoed  the  flow  fie  id  ana  offered  the  first  quantitative  data 
concerning  the  vortices  contained  within  the  field.  The  pressure 
measurement  of  the  field  produced  similar  quantitative  results 
concerning  the  position  of  the  primary  vortex  with  changes  in  aspect 
ratio.  These  measurements  also  yielded  information  concerning  the  flow 
field  pressure  distribution  with  changes  in  slot  aspect  ratio. 

The  boundary  layer  was  maintained  at  a  constant  thickness  by 
means  of  32  inches  (0.81  m)  of  number  12  rough  grain  sandpaper  and 
running  all  the  tests  at  15.75  ft/s  (4.8  m/s).  The  boundary  layer  profile 
was  measured  at  the  same  location  as  the  upstream  comer  of  the  slot. 
These  measurements  showed  a  shift  to  the  right  when  compared  to  the 
Soalding  flat  blate  boudary  layer  profile.  This  was  attributed  to  a  lack  of 
recovery  time  the  flow  had  between  the  sandDaper  and  the  location  of 
boundary  layer  profile  measurement. 


An  important  result  from  the  flow  visualization  was  documentation 
of  the  effect  of  reduction  of  the  slot  span.  All  indications  from  this 
investigation  show  no  significant  change  in  primary  vortex  location  or 
shape  when  the  span  was  reauced  from  38  inches  (0.97  m)  to  8  inches 
(0.20  m).  These  results  do  net  imply  that  three-dimensional  effects  don’t 
exist.  In  fact,  three-dimensionality  is  seen  in  some  photographs  where 
smoke  has  folded  over  on  itself.  Nevertheless,  the  time-averaged  results, 
and  in  particular  the  primary  vortex  location  seem  to  be  unaffected  by 
three-dimensionality.  The  photographs  of  the  flow  structure  with  various 
aspect  ratios  add  much  clanty  to  the  behavior  of  this  type  of  flow.  For 
example,  the  development  of  a  secondary  vortex  was  observed  to  occur 
at  aspect  ratios  of  1.75  and  2.0.  There  was  large  scatter  in  the  location  of 
the  primary  vortex  as  measured  from  the  flow  visualization  photographs. 
This  scatter  is  clearly  seen  by  the  error  bars  in  Figure  9.  The  scatter  is 
due  to  unsteadyness  of  the  flow  field,  and  possibly  also  to  three- 
dimensionality  of  the  flow.  Nevertheless,  the  trend  in  Figure  9  is  that  the 
vortex  center  moves  deeper  into  the  slot  as  aspect  ratio  increases  until 
about  AR  =  2.0,  beyond  which  the  location  of  the  primary  vortex  stays 
fairly  constant. 

The  pressure  measurements  inside  the  slot  offer  the  most 
enlightening  information  of  all  the  work  described  herein.  The  data 
shows  distinct  rings  of  constant  pressure  (isobars)  that  have  not  been 
measured  previously.  The  pressure  profile  at  aspect  ratio  0.5  shows  the 
vortex  center  to  be  very  close  to  the  traverse  limits  and  the  flow 


visualization  data  indicates  that  it  is  much  lower  and  out  of  range.  The 
lower  aspect  ratios  also  place  the  vortex  more  downstream  than  the 
higher  asoect  ratio  cases.  The  vortex  size  seems  to  scale  with  the  slot 
minimum  dimension,  but  this  does  not  explain  why  the  vortex  moves 
downstream  in  the  slot.  Regarding  the  vertical  location  of  the  vortex 
center  the  higher  aspect  ratio  cases  (above  AR  =  0.5)  maintain  the  same 
trend  as  the  flow  visualization  location  data.  Both  show  the  vortex  going 
deeper  into  the  slot  with  increasing  aspect  ratio.  They  both  show  the 
leveling  off  of  the  Drimary  vortex  vertical  location  at  an  aspect  ratio  of 
about  2.0. 

The  data  on  the  vortex  center  pressure  coefficient  is  more 
meaningful  when  discussed  in  terms  of  previous  work.  An  additional  plot 
was  made  which  superimposes  the  present  wind  tunnel  pressure 
coefficient  data  on  Ward's  water  tunne!  cavitation  index  data  (see  Figure 
20).  As  discussed  earlier,  if  the  observed  cavitation  by  Ward  was  purely 
vaporous  then 

o  ~  —  Cp  . 

This  implies  that  the  local  oressure  is  less  than  or  equal  to  the  vapor 
pressure.  F:gure  20  shows  that  an  order  of  magnitude  difference  exists 
between  the  two  sets  of  data.  This  was  the  same  discrepency 
commented  on  by  Holl  and  Billet  based  on  their  and  Ward's  cavitation 
index  data  ccmoareo  against  their  integrated  wind  tunnel  velocity  profile 
data  (Hoil  and  Billet  1936).  The  suggestion  is  that  the  cavitation 
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ccservec  m  me  vertex  center  s  cams  'y  gaseous:  w-  c"  ~c  es  mat me 
water  tunnel  sict  local  pressure  is  higher  than  vapor  pressure.  The 
external  pnctcgraens  succor;  this  conclusion.  Recai;  that  these 
pnctcgrapns  showed  that  mass  does  not  pass  into  the  primary  vortex 
cescite  the  mass  transfer  into  the  slot  through  the  entrance.  By  the  same 
argument  mass  can  not  leave  the  primary  vortex.  Hence,  if  gaseous 
buboles  aid  form  inside  the  vortex  by  gaseous  diffusion  they  would  be 
traccec.  Overall  the  present  data  confirms  Hoil  ana  Billet’s  suggestion 
that  the  cavitation  they  and  Ward  observed  was  gaseous. 

4. 1  Summary  of  Conclusions 

An  overview  of  the  experiments  conducted  for  the  slot  fluid 
mechanics  problem  is  in  order.  The  following  is  a  summary  of  the 
significant  conclusions  receded  here: 

1 .  The  turbulent  boundary  layer  at  the  model  entrance  location 
was  1 .2  times  the  slot  width.  The  profile  shape  was  similar  to 
Spalding's  profile  with  a  slight  rightward  shift  due  to  a  lack  of 


recovery  time  from  the  roughness  used  to  trip  the  boundary 
layer. 

2.  Flow  visualization  showed  no  change  in  the  two-dimensionality 
of  the  slot  fluid  mechanics  where  the  span  was  reduced  from 
38  inches  (0.97  m)  to  8  inches  (0.20  m). 


a . 

MM 


v,*.  a 
-  <•  "a*  v 

,-’^s’a’a“ 

a  "  *  **  M  *  O 

hv/.;V 
*.v.v>v 
•v-  -aV 


'  i 


•  *  i 


.  N  ^  N 


r.v.-'v'v 


V.*- V-  A  A 


B 


3.  r  ow  visualization  cemcnstratec  ‘"a:  "~e  s  c:  vortex  ccac.c" 
changes  with  ascect  ratio,  ana  sr.cwea  the  formation  of  a 
seccnc  vertex  for  asoect  ratios  of  2.3  ana  larger, 
i.  Mass  transfer  was  seen  to  occur  across  the  slot  entrance; 
nowever.  this  mass  transfer  c:a  not  extena  into  the  primary 
vortex  center. 

5.  Direct  pressure  measurement  showed  vortex  location  to  be  a 
function  of  asDect  ratio,  similar  to  the  benavior  snown  by  flow 
visualization. 


6.  Vortex  strength,  based  on  direct  pressure  measurement,  is 
shown  to  increase  with  decreasing  aspect  ratio.  This  and  the 
magnitude  of  the  pressure  coefficient  inside  the  vortex  lead  to 
the  conclusion  that  in  previously  documented  water  tunnel 
experiments,  the  ooserved  slot  cavitation  was  actually  gaseous 
rather  than  purely  vaporous. 

7.  Fiow  visualization  showed  that  the  direct  pressure 
measurement  did  not  alter  the  interna!  two-dimensional  fluid 


mechanics  of  the  slot. 
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Figure  A-3.  Cavity  Fiow  Visualization  at  AR  =  1.0 


Figure  C-4.  Pressure  Coefficient  Profile  An 


Figure  C-5.  Pressure  Coeflicient 


Pressure  Coefficient  Profile  AR  =  1.5 
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FUNCTION:  The  following  program  performs  several  tasks 

including : 

-  Recording  the  dynamic  head  of  the  wind  tunnel 
for  velocity  calculations. 

-  Recording  the  voltage  signal  from  the  hot  wire 
circuit. 

-  Controls  the  traverse  that  the  hot  wire  is 
mounted  on. 

FORTRAN:  The  following  program  was  written  for  use  in 
Microsoft  Fortran  version  3.2. 


EQUIPMENT:  The  following  equipment  and  software  was  used 

in  the  operation  of  this  program: 

-  Metrabyte  DASH-16  analog/digital  I/O 
expansion  board 

-  Metrabyte  24  bit  parallel  port 

-  Metrabyte  DASH-16  Fortran  Library 

-  MKS  high  accuracy  pressure  transducer  model 
398HD 

-  MKS  signal  conditioner  type  270B 

-  Hot  wire  and  bridge  circuit 

-  PSU  Mech.  Eng.  Wind  Tunnel 
X-Y  traverse  system 

OPERATION:  The  following  steps  are  needed  prior  to  run¬ 

ning  this  program: 

The  DASH-16  should  be  configured  as  follows: 

-  8  channel  differential  input 

(  WARNING  note  the  definition  of 
Metrabyte  'differential  input’.) 

-  unipolar  This  offers  a  0-5  vole  \ 

-  gain  set  at  2  J  range  of  measurement.  ! 

-  DMA  level  of  3 

-  base  address  set  at  768 


m 


>>>; 


pm 
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The  MKS  Signal  Conditioner  snouid  be  set  up 
as  f ol lows : 

-  sensor  range  set  at  1  torr 

-  zero  the  unit  ?.ocoro i r.g  to  manual 

-  initially  set  range  select  at  XI 

The  final  set  up  requirement  of  the  equipment 
is  to  have  the  pressure  and  hot  wire  signals 
placed  consecutively  in  that  order  on  the  A/D 
board.  This  is  for  later  convenience  and  is 
part  of  the  lab  procedure. 

Following  toe  completion  of  the  above  items 
and  noting  that  the  wind  tunnel  is  ready  for  a 
test  this  program  can  be  implemented. 

Just  prior  to  running  PROFILE  the  inlet  temp¬ 
eratures  of  the  wind  tunnel  and  the  atmospher¬ 
ic  pressure  should  be  noted  for  later  ref¬ 
erence.  These  will  be  needed  for  later 
calculation  of  the  wind  tunnel  velocity.  This 
is  the  time  to  place  the  hot  wire  probe  where 
you  want  to  start  the  boundary  layer 
measurement.  It  is  easeier  to  place  this 
origin  closest  to  the  surface  the  boundary 
layer  is  on  and  then  step  away  into  the 
freesream.  The  coordinates  of  the  pitot  probe 
relative  to  this  origin  are  also  required. 
PROFILE  will  place  the  hot  wire  at  that 
location  for  the  calibration  runs. 

Now  run  PROFILE.  Extra  care  should  be 
used  when  typing  in  the  step  size  and  limit 
of  the  measured  boundary  layer  profile.  For 
example  when  traveling  downward  in  a  boundary 
layer  both  the  y  step  size  and  the  y  limit 
should  be  negative  or  else  the  program  will 
crash.  Care  should  also  be  taken  so  that  the 
computer  will  not  run  the  hot  wire  into  any 
thing  solid.  Finally  the  user  must  record  the 
y  step  size  for  later  use  since  this  is  not 
done  by  computer. 


IMPLICIT  INTEGER *t2  (A-V) 

REAL  FREQ, DX. DY, SPEED 

DIMENSION  CALDA7 { 4 , 1C0C ) ,  DATA(4.5000) 

CHARACTERS  CALOME.  CAL  TWO,  DAT3T0R,  DATE,  TIME 

WRITE(*, ’ (A) ’ )  ' ***********CALIBRATION  INFORMATION**** 

WRITE( *, ’ ( A\ ) ' )  ’  ENTER  THE  FIRST  CALIBRATION  FILE  NAM 


WRITE;*, ’ (A\) ’ )  ’  ENTER  THE  SECOND  CALIBRATION  RILE  NA 

ME:  ' 

READ'  *, ’ (A) ' )  CAL TNG 

WRITE; *,  ’ ( A\) ’ )  ’  ENTER  CALIBRATION  SAMPLING  RATE/CHAN 
.  IN  HZ:  ' 

READ ( * , * )  SPEED 

WRITE; *, ’ (A) ' )  ’  ENTER  THE  FIRST  CHANNEL  TO  BE  SCANNED 
DURING ’ 

WRITE(*, ’ (A\) ’ )  ’  CALIBRATION:  ’ 

READ ( * , * )  STRTCAL 

WRITE; *, ’ (A) * )  ’  ENTER  THE  LAST  CHANNEL  TO  BE  SCANNED 
DURING ' 

WRITE; *,  ’ (A\) ’ )  ’  CALIBRATION:  ’ 

read;*,*)  fincal 

write; *, ' (A\) ’ )  ’  ENTER  THE  X  CALIBRATION  POSITION:  ’ 
READ ( * , * )  XCA— 

WRITE; *,  ’ (A\) ’ )  ’  ENTER  THE  Y  CALIBRATION  POSITION:  ’ 
READ;*,*)  YCAL 

WRITE;*, ’ (A\) ’ )  ’  ENTER  THE  NUMBER  OF  CALIBRATION  SETS 

.  * 

read;*,*)  calsets 

X0LD=0. 0 
YOLD=0. 0 


write;*, ’ ; A) ’ )  * *************sampling  information***** 
**********  ’ 

WRITE;*, ’ (A\)' )  ’  ENTER  THE  DATA  FILE  NAME:  ’ 

READ(*, ’ (A) * )  DATSTOR 

write;*, ’ ;a\) ’ )  ’  enter  data  sampling  rate/chan.  in  hz 

.  > 

READ;*,*)  FREQ 

WRITE;*, '(A)’)  '  ENTER  the  first  channel  to  be  scanned 
DURING’ 

WRITE;*, ' (A\) ’ )  '  SAMPLING:  ’ 

read;*,*)  sch 

write;*,  * ;a) ’ )  *  enter  the  last  channel  to  be  scanned 

DURING’ 

write;*,  ’  ;a\)  ’ )  ’  sampling-.  ’ 

READ ( * , * )  FCH 

WRITE;*, ’ (A\) ’ )  ’  ENTER  THE  Y  DIRECTION  STEP  SIZE:  ’ 
READ ( * , * )  DY 

WRITE;*, ’ (A\) ’ )  ’  ENTER  THE  MAXIMUM  Y  VALUE:  ’ 
READ;*,*)  YLIM 

WRITE; *, ’ (A\) ’ )  ’  ENTER  NUMEER  OF  SAMPLES  TAKEN  PER  PO 
SIT ION : ' 

READ(*,*)  NOS 
C 

C  CALCULATE  THE  NUMBER  OF  CONSECUTIVE  CHANNELS  TO  BE  SCANNED 
C  FOR  THE  CALIBRATION  AND  THE  SAMPLE  ROUNTINES 
C 

NUMCAL=NCCHAN ( STRTCAL , FINCAL ) 

N U M D  A T = NO CH AM ( SCH . FCH  ) 


SAMPLING:  ’ 


ENTER  NUMBER  OF  SAMPLES  TAKEN  PER  PO 
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C  THE  FIRST  CALIBRATION 


CALL  CALIBRATE  (  XCAL ,  YCAL ,  CALSETS  SPEED  ,  CALCNZ.  CTRTC 

.FINCAL,  '  . 

>  NUMCAL, CALDAT, XOLD,  YOLD) 

C 

C  THE  BOUNDARY  LAYER  MEASUREMENT  SUBROUTINE 
C 

CALL  SAMPLE  ( SCH, FCH, NOS, FREQ, DY, DATSTCR, NUMDAT,  YLIM, D 

AT  A, 

s  Aw j-< w  ,  i  vl.  j  y 


C  THE  SECOND  CALIBRATION 
C 

CALL  CALIBRATE  ( XCAL, YCAL,  CALSETS , SPEED, CALTWO, STRTCAL 
, FINCAL. 

>  NUMCAL , CALDAT ,  XOLD ,  YOLD ) 

STOP 

END 


INTEGERS  FUNCTION  NOCHAN ( START,  FINISH ) 


PURPOSE:  TO  CALCULATE  THE  NUMBER  OF  CONSECUTIVE  CHANNELS 
BETWEEN  THE  STARTING  AND  FINISHING  CHANNELS 


IMPLICIT  INTEGERS (A-V) 

IF  (START.LT. FINISH)  THEN 
NOCHAN= ( FINISH-START )+l 
ELSEIF  (START. GT. FINISH)  THEN 
NOCHAN=( 8- (START-FINISH) )+l 
ELSE 

NOCHAN=l 
END  IF 
RETURN 
END 


SUBROUTINE  CAL I BRATS (XCAL,  YCAL, CALSETS, SPEED, CALSTOR,  S 
TRTCAL , 

>  FINCAL, NUMCHN, CALDAT, XOLD, YOLD) 

IMPLICIT  INTEGER* 2  (A-V) 

REAL  SPEED, DISPLAY 
DIMENSION  CALDAT ( NUMCHN , 1000) 

CHAfiACTER*8  CALSTOR, RESPON 

WRITE ( *, '(/,/, A) ’ )  ’  ***#****^*****************#***** 


E  *’ 


WRITE( *, ’ (A) ’ ) 


*  BEGINNING  CALIBRATION  PRCCZDUR 


WRITE(*,  ’  (A,  /,/)’)  ’  *********************************** 
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MOVE  THE  HCT  **  I  EE  TO  THE  CALIBRATION  COORDINATES 
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DASHBS=768 
DMALZV-I 
INTLZ7=2 
E""IFLG  =0 

CALL  ADINIT( CASHES,  BMALEV, INTLEV, RTNFLG ) 

CALL  TRIGGER ( SPEED, NOMCHN) 

MODE=3+16 

SET  THE  NUMBER  OF  SAMPLES /'CHANNEL  TAKEN  PER  CALIBRATION 
POINT 

NUMPAS=1000 

RECLEN=NUMCHN*2CCC 

OPEN ( 10,  FILE=CALSTCR, ACCSS3=’ DIRECT’ , STATUS-’NEW’ , REC 
=RECLSN ) 

START  THE  CALIBRATION  SET 
DC  100  1=1, CALSETS 

WRITE  {  *,  '  {/,  A,  /}  '  )  ’  *********  **  ***  ******  ************ 
*******  ’ 

WRITE( *,  ’  ( A\ ) ’ )  ’  THE  CURRENT  CALIBRATION  SAMPLE  IS: 
WRITE(*, ’ (12, /) ’ )  I 

WRITE ( *,’(/, A\ )' )  ’  HIT  RETURN  WHEN  READY  TO  TAKE  CAL 
SAMPLES  ’ 

READ ( * , ’ (A) ’ )  RESPON 
TAKE  A/D  MEASUREMENTS 
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CALL  ADCCNV( MODE, STRTCAL. FINCAL, NUMPAS, CALDAT ( 1 , 1 ) , RT 


) I SPLAY  THE  FIRST  DATA  SAMPLE  OF  THIS  POINT 

WRITE ( 1 0 , REC  =  I )  ( { CALDAT f J , K ' , J  =  1 , NUMCHN )  ,  K-l ,  1000 ) 
WRITE  (  *  ,  ’  ’/,  A)  ’  ]  ’  ***  APPROXIMATE  DATA  VA.LL—S** 

CHAN  =  STRTCAL 

DO  30  K-l. NUMCHN 

DISPLAY--  CALDAT  ’  K.  1'  'CIO.  0) 

WRITE,'  *,  30}  CHAN,  DISPLAY 

FORMAT(/, IX. ’VALUE  ( VOLTS}  AT  CHANNEL  ’,11,’  =  ’ . F8 
IF  ( CHAN . LT .  THEN 
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o  o  o  non  o  o  o  o  o  *  *  one 


BOUNDARY  LAYER  ORIGIN. 


vc.  *  v,  r„*<  ^nn  nui  m  r.n. 


vnT  *  t 

AULJ-ALnu 

XNEW=0. 0 
YNEW=0. 0 

WRITE ( '*.  '{/,/,  A.  /}  ’  )  ’  ■**********£>.;£  OF  CALIBRATION ■»*■■« 

*:  '*  *  :fc  *  ;*  x  * 

WRITE (  *,  ’  {  /  ,  A)  ’  )  ’  ***RETUR1IIMG  TO  ORIGIN** 

? 

CALL  TRAVRS  ( KNEW ,  YNEW ,  XOLD ,  YCLD ) 

X0LD=0. 0 
YOLD-O. 0 
RETURN 
END 
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:N ,  YLIM, 

>  DATA, XOLD, YOLD) 

IMPLICIT  INTEGER *2  (A-V) 

REAL  FREQ, DX, DY, DISPLAY 
DIMENSION  DATA( NUMCHN , NUMPAS ) 

CHARACTER'S  DAT3TCR 
DX=0. 0 

RESET  THE  A/D  BOARD  TO  THE  SAMPLING  FREQUENCY 

cALi-,  TRIGGc,R(  FREQ ,  NUMCHN ) 

MODE-3 -13 
RTNFLG=0 

OPEN’  THE  SAMPLING  FILE 

REC=NUMCHN *2* MU MPAS 

OPEN ( 1 1 , FILE=DATSTCR, ACCESS^' DIRECT’ . STATUS = ’  NEW ’ 


*  -  VTTTT  rnfZ~7  r7  *  t 

.  i  .1 Vn.u,J 

OFF  ERROR. 


ALL-,  v^r  0.  1  -N  THE  Nc.XT  LINE  ALL^-ViATES  ROUND- 


vVSN 
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DELAY  LCCP 


)NTINUE 


C  TAKE  A/D  MEASUREMENTS  FGR  THIS  BOUNDARY  LAYER  POINT. 

C 

CALL  AD  CON  V(  MODE .  SCH ,  FCH ,  MUMP  AS ,  DATA(  1,1),  RTNFLO 
NR  ITE  (  11.  REC=I )  (  (  DATA  ( J ,  L )  ,  J  =  1 ,  MUMCHM  ;  ,  L  =  1 ,  MUMP  AS 

C  DISPLAY  THE  FIRST  DATA  SAMFLE  OF  THIS  POINT 


WRITE**, ’ (/, A) ’ )  ’ 


*** APPROXIMATE  DATA  VALUES** 


CHAM = SCH 

DO  90  K= 1 , NUMCHN 
D I SPLAY = ( DATA ( K .  1 ) /8 1 9 . 0 ) 

NR I TE ( * , 8  0 )  CHAN .  D I SPLAY 

FORMAT*/, IX, 'VALUE  (VOLTS)  AT  CHANNEL  ’,11,’  =  ’ ,  F8 

IF  (CHAN. LT. 7)  THEN 

CHAN -CHAN -1 

ELSE 

CHAN-0 

END  IF 

CONTINUE 

WRITE* *,  70) 

FORMAT  (  /,  IX,  ’************************************** 


(*'  ,  /; 


0  CHECK  PROGRESS  IN  MOVING  THOURGH  THE  BOUNDARY  LAYER 
C 

Ir  (  I .  i3r.  I  STEPS)  GOTO  ICO 


X'  I E  W = XCLD + D  X 


J-  : 


RETURN 

END 
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WILLIAM  G.  BLAIR 
10/17/36 


C  PURPOSE:  To  convert  binary  files  with  integer*2  data  and 
C  write  the  data  to  a  second  formated  file. 

C 

C 

C  FORTRAN:  The  following  program  was  written  for  use  in 


Microsoft  Fortran  version  3.2. 


FILI 


IMPLICIT  INTEGERS  (A-Z) 

CHARACTER* 1 4  BFILE, DFILE 
DIMENSION  DATA( 4, 1000 ) 

WRITE ( * , ’ ( A\ ) ’ )  *  ENTER  THE  NAME  OF  UNFORMATED  BINARY 

- >  » 

READ( *, ’ (A) ’ )  BFILE 

WRITE (  *,  * ( A\ ) ’  )  *  ENTER  THE  NAME  OF  FORMATED  DECIMAL  F 

ILE - >  ' 

READ(*, ' (A) ’ )  DFILE 

WRITE( *, ’ ( A\) ’ )  ’  ENTER  THE  NUMBER  OF  RECORDS  IN  FILES 


■> 


EL’ 


READ ( * , * )  NREC 
WRITS ( *, ’ (A\) ' J 

READ ( * , * )  CHAN 
WRITE( *, ' (A) ’ ) 


ENTER  THE  NUMBER  OF  CHANNELS- 


WRITE ( *, ’ ( A\ ) ’ )  ’ 
- >  ’ 

READ( *, *)  PNTS 
LENGTH=CHAN*PNTS*2 


ENTER  THE  NUMBER  OF  POINTS  PER  CHANN 
(IF  POINTS  .GT.  1000  -  ALT  PROGRAM) 


ENGTH 


WRITE ( *, 

J  (/,/ 

.  A)  ’  ) 

WRITE ( *, 

’  (A)  * 

)  »  * 

WRITE ( *, 

’  (A)  * 

)  *  * 

WRITE ( *, 

’  (A)  * 

)  *  * 

WRITE( *, 

’  (A)  * 

)  '  * 

WRITE (*, 

’  (A)  ’ 

)  ’  ’ 

WRITE ( *, 

90) 

OPEN ( 10, 

FILE= 

BFILE 

OPEN (20, 

V- 

A.  A_<  - 

DFILE 

DO  10  K= 

1 ,  'iRE 

C 

***************** 


THE  FIRST  POINT  OF  EACH 
CHANNEL  FROM  EACH 
RECORD 


*' 
*  ’ 
* ' 


*****************' 


10 


READ( 10, REC-X)  ( (DATA( I, J) , 1=1, CHAN) , J=l, PNTS) 
WRITE ( * , ICO)  (DATA(I, 1), 1=1, CHAN) 

DO  10  J=l, PNTS 

WRITE ( 20, 100)  (DATA(I, J), 1=1, CHAN) 

CONTINUE 

WRITE( *, ’ ( A) ’ )  ’  ******  EOF  REACHED  ******* 
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io  calculate  one  least  square  curve  :::s  of  the 
ca i i o r at i on  cata  or  rrccuce  a  rev  ifile  >■  — 

of  the  boundary  layer  profile  points. 

The  following  rrcgram  performs  several  tasks 
including : 

-  Reads  and  averages  each  record  of  the  given 
file.  The  new  mean  value(s)  is  either  a 

calibration  velocity  pair  or  a  boundary  layer 
point . 


For  calibration  data  the  mean  values  of  the 
dynamic  head  in  the  first  column  of  the  array 
are  used  to  calculate  the  wind  tunnel 
velocity  for  that  point. 

For  calibration  data  a  least-squares  curve  fi 
of  order  5  is  performed  on  the  velocity  vs 
voltage  data.  Both  the  polynomial  coeffi¬ 
cients  and  the  velocity  vs  voltage  data  are 
stored  in  a  new  file. 

For  boundary  layer  data  the  mean  hot  wire 
voltages  are  mapped  onto  their  velocity 
counterparts  with  the  calibration  data  curves 
The  calibration  curve  used  is  a  linearly  time 
weighted  combination  of  the  first  and  second 
calibration  sets. 

For  boundary  layer  data  the  position  vs 
velocity  curve  is  stored  in  a  new  file. 


The  following  program  was  written  for  use  in 
Microsoft  Fortran  version  3.2. 


This  program  dees  not  do  ail  of  its  functions 
in  one  session.  The  calibration  files  need  t 
be  curve  fit  prior  to  processing  the  boundary 
layer  points.  This  requires  SLAYER  to  be 
run  three  separate  times  for  each  boundary 


The  actual  operation  of  this  nrogram  is  self- 
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C  are  ir  rneti e r s / s ec  -irici  vc  ltd.  T^.-^  3  _  ; 

C  1  eye  r  d e~  e  3  ^ a  i r. c. r. ^ s  iii .1  rr.e  ~  r  r  z  j  *■?  i- 

C  English  and  Metric  ur.:^  ver-i 

C  for  convenience  since  cur  wind  tunnel 

C  velocities  are  offer,  found  ir.  r.  s  ant  t.- 

C  traverse  used  for  hot  wire  movement  has 

C  of  the  data  is  done  with  other  software. 

C 

IMPLICIT  INTEGER*2  (A-H.O-Z) 

INTEGER  WAVE 

HEAL  MEAN,  VAR,  A,  ANS,  X,  V.  PP.CFIL,  HHO,  CONST.  rP.ns, 

T, DY, YY 

<p  r  ■>  ro 

CHARACTER*20  FILS1, FILE2, FILES, FILE4 
DIMENSION  DATA( 1000, 3 ) , MEAN ( 3 ) , PROFILE  50, 3), A( i: 

(100) 

3  ,  Ad  S  (  2  0  )  ,  Cl  (  10)  ,  02  (  i.0) 

COMMON  X(100), Y( 100) 


FIT=0 

WRITE (A\) ’ )  ’  ENTER  THE  NAME  OF  THE  DECIMAL  DATA  F 


00 ) , VAR 


READ ( *, ’ ( A) ’ )  FILE1 

WRITE ( *, ’ ( A\) ’ )  '  ENTER  THE  NUMBER  OF  CHANNELS  IN  THE 

•v  > 

- ; 

READ  (  * ,  * )  CHAN 

WRITE( *, ' (A\) ’ )  ’  ENTER  THE  NUMBER  OF  POINTS  PER  CHANN 
- >  ' 

REAP ( * , * )  PNTS 

WRITE( *, ' (A\) ’ )  ’  ENTER  THE  NUMBER  OF  RECORDS  IN  THE  F 
- >  * 

READ (  * ,  *  )  WAVE 

WRITE; *, ’ (A) ’ )  ’  DO  YOU  WANT  LEAST  SOU ARE  FIT  OF  THE  1 
iANNEL  ’ 

WRITE! ’ (A\) ’ )  ’  vs  2nd  CHANNEL  (1/0)  - 


READ!*.*)  FIT 


TC  '  r*  -  x,  «  /  a  »  >  >  t  *7 1 


OPEN( 11, FILE=FILE1. STATUS- 'OLD’ ) 
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write: *, 

(  A  )  ’  y 

•  ENTER  THE  NEW  DATA  FILE 

NAME - 
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f'ijX  L  *  H  [  *■  , 
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C  PHCFIL  IS  AN  AVERAGE  VOLTAGE  BASED  ON  AN  A/D  RANGE  OF  0-5 
C  VOLTS  USING  THE  12  BIT  BOARD. 

r* 

??.CF  IF'  I .  J  '  =ME A- 1 '  J  ;  *'5.0/4035.0 


CT>  T  r  nr  or  Tirr  T7  AT  7T7  ^  a  ■*“*  T  Mnn  /"'XT  *n* 

mu  nvnnriun  /m-.onn  xnn  O'^rcnnd 

c 

DO  17  1=1. WAVE 

WRITE  ( •* ,  S3  )  I ,  (  J ,  J  =  1 , 4 ) 

WRITE ( *, 9S )  ( PROF IL ( I , J ) , J= 1 , CHAN ) 

i  *->  Tvt:tT7 

C 

C  CHECK  TO  SEE  IF  A  LEAST-SQUARE  CURVE  FIT  IS  DESIRED 
Ir  (  j  =T .  =  .  1 )  TH=.N 

r'*r--.nv  ;  ia\/  ;  n-i-nn.  lor.rnJ  i  ijwv^riL,  rrlnnnunn  (  crnng 

) - >* 

READ  (*,*)  PRES 

WRITE ( *, ’ ( A\) ’ )  ’  ENTER  THE  CURRENT  ROOM  TEMPERATURE 
(deg  C)  — >' 

READ ( * , * )  TEMP 
C 

C  CONST  CONVERTS  VOLTS  TO  PRESSURE  UNITS 
C  CONST  =  ( Pa/mmHg ) *( mmKg/Volts ) 


C  RHO  IS  DENSITY  kg/m**3 
C  GAS  CONSTANT  FOR  AIR  E=23S.3  J/(kg*K) 

RHO  = .'  PRES  *  1 0 .  0*C0MST)  /(  (  TEMP -*-273  .  15}  *235. 3) 


WRITE  (*.■*;•  ' . 


.ST  CHANNEL  CF  PRESSURE  IS  CONVERTED  TO  VELOCITY 


PROFILd,  1 )  =SQRT  (  2  .  0*PRCFIL(  I,  1 )  *( CONST/ 10.  0) /RHO 
CONTINUE 
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INTEGER*2  (A-H.O-Z) 


[MFLIGIT 
IEAL  MEAN 

^  7  \4T'.J  C  7“  Q \J 


sA  .  rt  \  _  L/C  .. 


MT 


DO  10  1=1, CHAN 
MEAN ( I ) =0 . 0 

T'O  s.  "111'C 


MEAN ( I ) =MEAN ( I ) +FL0 AT ( DATA ( J ,  I ) ) 
CONTINUE 

MEAN ( I ) =MEAN ( I ) /FLOAT ( PNTS ) 
CONTINUE 
RETURN 
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-URPCSR:  TO  PER? 


CALCULATE  THE  MOMENTUM  TI- 

V'T’'  *  r  r  \ 

O  4.  J  Lftit  :  Lft  :  Z.i* . 


ORATION  IN  ORI 

»  i_  ^  *0  l* r*  4^  ^  j  , 


me  following  program  was  written 
m;  ^r-'^cf f  "orr’”'  o v_'  3  n 


;e  m 


OPERATION:  The  only  needed  information  on  the  boundary 

layer  is  the  freestream  velocity  and  the  number 
of  points  in  the  file.  Otherwise,  this  program 
is  self-explanatory . 


IMPLICIT  REAL*8 ( A-H ,  0-Z ) 

INTEGERS  PNTS 
REAL *8  NU 
CHARACTER* 1 4  FILE1 
DIMENSION  U(100), Y(100) 

DATA  NU/1 . 508D-05/ 

WRITE ( *, ’ ( A\) ’ )  ’  ENTER  THE  MEAN  VELOCITY  (m/s) 


■>  ’ 


READ ( * , * )  UFREE 

WRITE (  *,  ’  ( A\ )  ’  )  ’  ENTER  THE  BOUNDARY  LAYER  FILE  NAME— 


RSAD( *, ’ (A) ’ ) FILE1 

WRITE ( *,  ’ ( A\) ’ )  ’  ENTER  THE  NUMBER  OF  POINTS  IN  THE  FI 


•>  ’ 


READ ( *, * )  FNTS 

OPEN( 10, FILE=FILE1, STATUS=’OLD’ ) 
DO  10  1=1, FNTS 
READ (10,*)  Y ( I ) , U ( I ) 


ALL  Y’S  HAVE  0.1  INCH  ADDED  TO  SHIFT  DATA  TO  THE  ACTUAL 
SUr„“  ACu  THn  BOUNDAn. L  LAYnn.  IS  FOLLOWING.  THEN 
CONVERT  Y’S  TO  METERS  FROM  INCHES. 


Y ( I )  =  ' Y ( I ) +  0 . 1 } * 0 .  0254 
CONTINUE 

i  iC) 


AREA  UNDER  CURVE  IS  IN  IT 1 1 


:  Ls  -if". by 


SUM=0. DO 

CALL  INTEGRAL  FUNCTION  FOR  FIRST  POINT 


yyy:/. 
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velcc:;y  13  required.  This  is  found  by  run¬ 
ning  MCMEN  specifically  made  for  that  purpose. 
Gnce  this  information  is  avalible  then  the 
operation  of  NCNDIM  is  sel f -exp iainatory . 


IMPLICIT  REAL *8  (A-H.O-Z) 
REALMS  NU 

CHARACTER *<20  FILEl ,  FILE2 
DATA  NU/1. 508E-05/ 


V  IS  THE  WALL -FRICTION  VELOCITY  PREVIOUSLY  CALCULATE! 


WRITE ( *, ' (A\) ’ )  ’  ENTER  V* 


READ (  * ,  * )  V 

WRITS( *, ’ (A\ ) ’ )  '  ENTER  THE  BOUNDARY  LAYER  FILE  NAME-- 


READ ( * , ’ ( A) * )  FILEl 

WRITE ( * , ’ ( A\ ) ’ )  ’  ENTER  THE  NUMBER  OF  POINTS  IN  THE  FI 


READ ( * , * )  NUM 

WRITE ( * , ' ( A\ ) ' }  ’  ENTER  THE  NONDIMENS ICN ALICES  SLAYER 


READ' *, ’ ( A; ’ :  FILES 

OPEN ( 10, FILE=FILE I , STATUS- ’ OLD ' } 


10  10  1=1, NUM 
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IKE  LOCAL  PRESSURE  IN  A  2-D  SLOT  0?  VARYING 
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-  Recording  the  dynamic  head  of  the  wind  tunnel 
for  velocity  calculations  and  ncnd imens icnai - 
is an  ion  of  the  pressure  difference. 

-  Recording  the  pressure  difference  between  the 
pressure  at  the  center  of  slot  base  and  the 
static  pressure  of  the  free  stream. 

-  Records  the  pressure  difference  between  the 
static  pressure  of  the  free  stream  and  the 
pressure  of  the  slot  probe. 

-  Calculates  the  absolute  X  and  Y  coordinates 
relative  to  the  up  stream  corner  at  the  slot 
entrance  and  records  the  calculated  location 
with  each  pressure  data  set. 

-  Calculates  the  relative  X  and  Y  coordinates 
for  interactive  use  between  the  user  and 
‘Dsriiiincil .  These  coordinates  t^ell  ^he  user  ishe 
X  and  Y  values  for  the  traverses;  since  they 
need  to  be  set  by  hand. 


rCRTRAN :  The  following  program  was  written  for  use  in 
Microsoft  Fortran  version  3.2. 


The  following  equipment  and  software  was  used 
in  the  operation  of  this  program: 


c 

-  Metrabyte  DASH-16 

an  alog/ digital 

1/0 

expansion  bears 

-  Metrabyte  24  bit 

parallel  port 

-  Metrabyte  DASH- 13 

~ y  0  y>  I ,  ’  h  r  R.  r 

V 

-  MRS  high  aoouraov 

d  ;  *  *  e  r  e  n  t  i  a  1  t 

ressure 

^ r^r.od’jc - r  mod*?! 

333HD 
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-  MRS  signal  tend  it 

1 0  r.  e  r  type  2  ?  C  3 

’  I  pvt  -  >J 

The  fallowing  step 

rj  0  ''a  .0 ^  v. 

'  ^  *•  ^ 1 4 

nmg  this  program; 

c 

-  The  DASH-16  should 

be  configured 

as  follows 
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*c vs : 

-  sensor  range  ser  ar  i 

-  zero  the  unit  acccrdir 

-  initially  sen  range  s* 
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Configure  the  model  accordingly  f< 


tne  new 


test.  This  includes  the  following: 

-  mounting  the  correct  side  walls  for  the 
test 

-  centering  the  sloe  pressure  probe  in  the 
middle  of  the  2-D  grid  to  be  scaned  (use 
calipers ) 

-  record  the  absolute  and  relative  coordin¬ 
ates  from  the  model  and  the  traverses 
respectively 

Following  the  completion  of  the  above  items 
and  noting  that  the  model  is  correctly  config¬ 
ured  for  a  test  this  program  can  be 
imp  1 e men ted . 
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vIlis'C  nr*'  o'*'  oo  nv-inni^g  ?Hn3sL*r*0  tins  iniatw 
eratures  of  the  wind  tunnel  and  the  atmespner- 
ic  pressure  should  be  noted  for  later  ref¬ 
erence.  These  will  be  needed  for  later 
calculation  of  the  wind  tunnel  velocity. 

The  use  of  the  relative  and  absolute 
coordinate  systems  needs  some  e>rplanaticn. 

The  absolute  coordinates  are  the  recorded 
positions  of  the  slot  pressure  probe  with 
reference  to  the  up  stream  corner  at  the 
slot  openning. 
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dative  coordinates  are  only 


)  r  user 


benefit.  They  refer  to  the  two  manual  traver¬ 
ses  and  their  coordinates.  Note  that  the  use: 
must  input  the  X  and  Y  for  both  of  the  two 
traverses  since  they  will  not  necessarily  be 
the  same.  nice  tne  inputs  nave  teen  typed 


:raverses 


i  rther  than  two  hit  the  return  hey  f cr 


:  n c  e  w  n  e  ^reverses  a  r 1 


set.  However,  the  user  will  be  busy  connect¬ 
ing  and  disconnecting  the  tubing  on  the  pitot 
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CHARACTER*20  DAT3TGR 

COMMON  XRACZ ,  XCEFF,  XCSNN ,  '/FACE,  YCENN,  YCI 


WRITE ( * ,  ’(A)’)  ’ *************SAMPLING  I M FORMAT I C  N  *  *  *  *  * 

X  :k  :<  *  *  'X  x  *  * 

WRITE  (  * ,  ’  (  A\ }  ’  )  ’  ENTER  THE  BATA  RILE  NAME :  ’ 

READ (  k ,  ’  (A)  ’  )  DA.TSTCR 

write  (A\)  ’  )  ’  ENTER  DATA  SAMPLING  RATE /CHAN .  IN  HZ 

> 

READ (  * ,  * }  FRED 

WRITS (*,  '  (A)  '  )  ’  ENTER  THE  CHANNEL  TO  3E  SCANNED  DUF.IN 

WRITE ( *, ’ ( A\ ) ’ )  ’  SAMPLING :  ’ 

READ ( * , * )  SCH 

WRITE (  *, ’ { A\ ) ’ )  ’  ENTER  THE  NUMBER  CR  GRID  POINTS:  ’ 
READ ( *, *)  NPNTS 

WRITE (  *, ’ ( A\ ) ’ )  ’  ENTER  THE  ABSOLUTE  RACE  CENTER  X  (in 

\  .  > 

READ  (■*,*)  XRACE 

WRITE  (  *,  ’  ;  A\ }  ’  )  ’  ENTER  THE  ABSCL'JTE  RACE  CENTER  Y  (in 


READ ( *,*}  YRACE 
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READ  XCE'NN 
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WRITE (*,’ (A\ )’ )  ’  ENTER  THE  RELATIVE  FAR  CENTER  Y  (in) 
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SUBROUTINE  SAMPLE  {  Nr  MTS ,  SCH ,  FCH ,  MUMP  AS  ,  FRE< 
IMPLICIT  INTEGER  *2  (A-V) 

iviiiiL  P  ml*  <j/.  j  U  1  or'  i  f  a  i  *. 

DIMENSION  DATA ( NUMCHN , MUMP AS ) 

CHARACTER* 8  DATSTCR, JUNK 
COMMON  XFACS , XCEFF , XCSNN , YFACZ , YCENN , YCEPF 


'•rry  m  r\  '/■  __  4  /  4  £  *  HP.  V—  5  /  *  £  CD  rnU  'T»0'  TT  ^  TD  C  ^ 
*-*^*ur  A-’i/ib  ru'i t  —  -1/-.0  r di ic  io  r  j.rt~  *. 

Y  9  kook. 

C  CALCULATE  GRID  SIZE  (i.e.  7x7  or  14x14) 

NUMX= .NT  (  SoF.T  {  -  LCAT  ( NFNT3 )  )  ) 

C 

C  SET  STEP  SIZE  AND  JUMP  TO  START  X  AND  Y 
C 

XD=0. 0625*7. 0/ ( FLOAT (NUMX) ) 
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Y = YF ACE -X JUMP 
XCEN =XCENN - X JUMP 
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*  —  /  i.  a—#  +  ?  Au  j. 

YCSN = YCENN - XJUMF 

Y  CEF = Y CEFF-XJUMP 
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WRITE( *, 80)  CHAN,  DISPLAY 
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IP  (CHAN.LT.7)  THEM 
CHAN=CHAN+1 
ELSE 
CHAN=0 
END  IF 

SO  CONTINUE 

WRITE ( *, 70) 

202  CONTINUE 
C 

C  RESET  THE  X  COORDINATES  FOR  A  PASS  ON  A  NEW  ROW 
C 

T =XF  ACE-XOU^^ 

XCE  N = XCZ N N - X J  U MP 
XCEF =XCEFF +XJUMP 
201  CONTINUE 
WRITS ( *, 70) 

70  FORMAT (  /  IX  ’  ^*^*************j!**x**:**^*±***'i****;**-Jc*xl 

,  / ) 

CLOSE ( 11) 

RETURN 

END 
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following  program  was  written  for  use  in 
csoft  Fortran  version  3.2. 
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